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P RE FACE 


Under  the  sponsorship  of  the  Tactical  Technology  Office  of  the 
Defense  Advanced  Research  Projects  Agency,  Rand  has  been  engaged  in 
the  development  of  hydrodynamic  design  criteria  employing  concepts  of 
boundary- layer  control. 

In  the  absence  of  a comprehensive  theory  for  the  location  and 
properties  of  the  transition  region  between  laminar  and  turbulent  flow, 
much  of  the  analysis  and  design  is  empirically  based.  However,  a 
secondary,  but  still  important  issue  can  nevertheless  be  approached 
from  an  analytic  viewpoint:  the  effect  of  boundary-layer  control  on 

the  separation  of  a laminar  flow.  Using  simple  analytic  methods  and 
results  from  boundary-layer  theory,  this  report  compares  the  effective- 
ness of  different  types  of  boundary-layer  control  in  delaying  laminar 
separation. 

The  report  should  be  useful  to  hydrodynamicists , designers  of 
submersibles , and  others  engaged  in  applying  fluid  mechanical  method- 
ology to  the  improvement  of  underwater  vehicle  performance. 


SUMMARY 


Boundary-layer  control  provider  a potentially  effective  method 
for  reducing  the  .kin  friction  of  submersible  vehicle..  Once  the  flo* 
1.  lamlnarlzed,  the  focus  often  ahlfta  to  the  problem  of  delaying  the 
oneet  of  lamlm.r  aep.r.tlon,  .her.  the  aurf.c.  .hearing  atr...  vanl.he. 
and  the  boundary  layer  1.  liberated  from  the  at.blllaln,  effect  of  the 
•all.  Tran. 1 cion  to  turbulent  bound.ry-l.yer  flo.  often  occur,  a.  . 

nsequence  of  this  separation,  because  of  the  enormous  rate  of  growth 
In  disturbances  In  the  separated  boundary  layer. 

The  effect,  of  suction  on  the  delay  of  separation  are  well  known- 
ana lyses,  computation.,  and  „...„r.„e„t,  ov.r  th.  p,.t  ,e„nty  ye„, 

have  ahown  that  laminar  aep.r.tlon  can  be  delayed  almost  Indefinitely, 
if  the  additional  complealtle.  and  coat,  of  auction  are  tolerable. 
However,  th.  effects  of  aurf.c.  ha.tlng  on  delay  of  aepar.tlon  ere 
leas  well  known— there  are  no  measurements  available,  end  there  ere 
no  directly  applicable  analyses  or  computations. 

To  fill  thia  gap,  wa  pre.ent  an  analysis  of  the  minimum  surface 
overheat,  Tw  - that  will  delay  separation  for  a prescribed  adverse 
pressure  gradient  in  water.  Laminar  boundary  layer,  usually  separate 
when  the  pressure  begins  to  rise  aft  of  the  maximum  diameter  on  a body 
of  revolution.  The  adverse  pressure  gradient  (dp/dx  > 0)  that  can  be 
sustained  without  separation  is  small  in  the  absence  of  heating  or 
suction.  However,  slight  alterations  of  the  velocity  profile  near  the 
wall  can  have  strong  effects  on  separation. 

The  analysis  is  for  a classic  type  of  boundary-layer  flow,  the 
Falkner-Skan  wedge  flow  corresponding  to  negative  valuer  6.  We 
perform  an  asymptotic  analysis  of  such  flows,  including  the  effects  of 
viscosity  variation  with  temperature.  The  energy  and  momentum  equa- 
tions are  coupled  through  this  variation,  and  we  employ  a high  Prandtl 
number  approximation  to  develop  systematic  approximations  to  the  equa- 
tions of  motion  and  energy.  We  find,  as  might  be  expected  intuitively, 
that  the  effect  of  heating  and  viscosity  variation  is  localized  to  a 


layer.  The  primary  effect  of  this  layer  on  separation  is  to  provide 
a "slip"  velocity  for  the  outer  main  parts  of  the  flow.  This  slip 
velocity  enables  the  ouLer,  shear-layer-like  part  of  the  flow  to  sus- 
tain a more  adverse  pressure  gradient  than  it  could  in  the  absence  of 
heating. 

However,  although  heating  is  shown  to  delay  separation,  the  mag- 
nitude of  its  effect  is  found  to  be  small,  particularly  if  practical 
values  of  wall  overheat  are  considered.  We  show  that  the  magnitude 
of  the  acceptable  adverse  pressure  gradient  parameter  increases  by 
about  25  percent  for  a 40°F  overheat.  For  a wedge  flow,  this  means 
that  U ~ x-0*11  rather  than  U ~ x 0,0^°  as  occurs  when  there  is  no 
heating. 

Comparisons  with  available  results  tor  delaying  separation  using 
suction  show  that  a suction  velocity  ratio  of  less  than  0.0001  would 
have  a comparable  effect  in  maintaining  an  attached  flow. 

Our  results  suggest  that  the  effects  of  surface  heating  on  the 
delay  of  separation  are  probably  real  and  significant,  but  heating 
does  not  have  the  almost  unlimited  capacity  for  maintaining  an  attached 
flow  that  suction  appears  to  possess.  In  addition,  it  should  be  empha- 
sized that  neglecting  the  effects  of  heat  transfer  in  prediction  of  the 
location  of  laminar  separation  produces  design  procedures  that  are 
essentially  conservative. 
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I. INTRODUCTION 


We  ere  interested  in  the  feasibility  of  delaying  the  separation 
of  a laminar  boundary  layer  that  is  subject  to  an  adverse  pressure 
gradient  in  water.  The  significance  of  maintaining  an  attached  boundary 
layer  in  regions  of  adverse  pressure  gradient  is  that  transition  to 
turbulent  boundary-layer  flow  often  occurs  by  laminar  separation  and 
turbulent  reattachment  slightly  aft  of  the  point  of  maximum  diameter, 
and,  as  body  size  increases,  the  relative  portion  of  body  surface 
under  turbulent  flow  Increases.  It  is  then  desirable  to  consider 
additional  means  for  controlling  the  location  of  laminar  separation. 

Two  such  methods  for  delaying  separation  in  the  presence  of  an 

adverse  pressure  gradient  are  suction  and  wall  heating.  The  effects 

of  such  are  vrll  known.  Schlichting^^  cites  measurements  made  by 

Prandtl  over  seventy  years  ago  that  verified  the  powerful  effect  of 

suction  on  delaying  separation.  The  effects  of  wall  heating  are  not 

so  well  known  and  do  not  yet  appear  to  have  been  measured.  However, 

(2) 

the  calculations  of  Wazzan  et  al.  for  a laminar  water  boundary 
layer  of  the  Falkner-Skan  type  subject  to  an  adverse  pressure  gradient 
indicate  that  surface  heating  in  water  could  maintain  an  attached 
boundary  layer,  which  would  have  separated  in  the  absence  of  any  tem- 
perature difference  between  surface  and  free  stream.  Uazzan's  calcula- 
tions for  a Falkner-Skan  S « -0.198,  which  is  the  maximum  negative 
value  for  6 for  unseparated  flow,  also  indicated  that  moderate  surface 
heating  in  water  (Ty  - T^  < 40T,  and  T^  - 60*F)  does  not  increase  the 
minimum  critical  Reynolds  number  as  much  as  it  does  for  flat-plate 
flow.  This  result,  although  not  hearing  on  the  subject  of  separation, 
strongly  suggests  that  even  if  surface  heating  could  be  used  to  main- 
tain an  unseparated  flow,  the  stability  characteristics  of  this  flow 
would  still  make  it  susceptible  to  transition  near  the  pressure  minimum. 
The  situation  is  probably  different  with  suction,  where  experience 
suggests  that  the  delay  of  separation  also  leads  to  the  delay  of 
transition. 
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In  this  report,  we  analyze  the  effect  of  wall  heating  on  the  de- 
lay of  separation  for  the  simplest  important  boundary-layer  flows  with 
adverse  pressure  gradient:  the  Falkner-Skan  flows  with  negative  B and 
negative  m,  where 


and 


Using  the  idea  of  local  similarity,  it  is  then  possible  to  extend  the 
results  developed  here  to  a more  general  boundary-layer  flow  with 
reasonable  accur.cy. 

We  are  primarily  interested  in  sharper  estimates  of  the  effect 
of  heating  in  order  to  compare  tire  efficacy  of  suction  versus  heating 
for  the  prevention  of  separation.  For  that  comparison,  analyses  based 
on  the  Falkner-Skan  flows  are  useful  because  the  effects  of  suction 
for  such  flows  have  already  been  examined  carefully  by  Terrill ^ and 
Nickel . ^ 

The  question  we  consider  is  the  following:  How  much  adverse  pres- 

sure gradient  will  a heated  laminar  boundary  layer  in  water  sustain 
without  separation?  In  the  language  of  Falkner-Skan  flows,  the  ques- 
tion than  becomes:  What  is  the  largest  negative  value  of  B for  an 
attached  flow  that  a given  wall  overheat  AT  = T^  - Toj  will  permit? 

The  analysis  is  a systematic  approximation  to  the  boundary-layer  equa- 
tions, based  on  the  properties  of  high  Prandtl  number  solutions,^ 
and  the  existence  of  a thermal  boundary  layer  that  is  considerably 
thinner  than  the  velocity  boundary  layer. 

The  Prandtl  number  of  water  at  608F  is  about  8 and  all  the  exist- 
ing evidence  is  that  this  is  large  enough  for  the  validity  of  methods 
that  are  based  on  high  Prandtl  number  arguments. 

For  water  flows  in  the  range  of  40°F  to  100°F,  the  principal  de- 
parture from  constant  property  flow  is  due  to  viscosity  variation. 

The  density,  thermal  conductivity,  and  specific  heat  do  not  vary  apprec- 
iably in  this  range,  and  the  most  important  phenomena,  because  of  the 
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gradient  of  viscosity,  are  quite  well  represented  by  this  model  in 


which  p,  k,  and  cp  are  constant,  and  the  variation  of  viscosity  with 
temperature  is  preserved.  There  exists  an  analogy  between  wall  heat- 


ing and  suction  for  a variable  viscosity  fluid.  For  example,  the 
compatibility  condition  at  a wall,  for  both  suction  and  wall  heating, 


w 3y 


-£♦ 


n.2  3p  ~ 

3 u , w 3u 
w 3y2  3y  oy 


or,  since  p = p(T) 5 


_ / d^J  3T 
~ \ dT  3y 


\ 

/wj 


dp  . 3 u 

- j + p — r 

dx  w » 2 

3y 


For  water,  dp/dT  < 0;  for  suction,  < 0.  Therefore,  since  the  tem- 
perature gradient  dT/dy  is  negative  if  Tw  > Tm,  the  effect  of  surface 
heating  would  be  qualitatively  similar  to  that  of  suction.  The  sim- 
ilarity is  only  qualitative — suction  has  a more  powerful  and  global 
effect  on  a boundary  layer  because  of  its  impact  on  boundary-layer 
profiles  and  the  growth  of  the  boundary  layer.  The  effect  of  heating 
is  localized  to  a smaller  region  (which  nevertheless  is  an  important 
one  for  stability),  and  there  is  little  effect  on  boundary-layer  growth. 


I 


W 
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II.  ANALYSIS 

LARGE  PRANDTL  NUMBER  APPROXIMATIONS 

The  boundary- layer  equations  are  written 


Momentum: 


„ Du  Du 

P u 7T-  + p v r— 
“ Dx  ' 1X1  Dy 


D_ 

Dy 


("«>  I 


0.) 


Energy: 


DT 

c u 

° p Dx 


, _ dt  . d2t 

+ P-V  37  - K*  7T 
Dy 


Continuity : 


D__ 

Dx 


(ur^ ) 


+ (vrj)  = 0 


(2) 

(3) 


where  j = 0 corresponds  to  plane  flow, 

j = 1 corresponds  to  axial-symmetry. 

If  the  new  variables 


e = 


dx 


and 


n 


y 


are  used  to  replace  x and  y,  and  if 


g = 


and 


Df  _ u 

Dn  = u 

e 


are  chosen  as  dependent  variables,  then  the  coupled  equations  of  vari- 
able viscosity  ?alkner-Skan  flow  are  obtained  provided  that 


where  N = Vi/MOT,  Pr^  * M^Cp^/k^,  and  primes  denote  differentiation  with 
respect  to  n*  ihese  equations  must  be  supplemented  by  the  appropriate 
boundary  conditions. 


f(0)  - f'(0)  - 0 , 

g(0) 

f'(“)  - 1 , 

g(°°) 

The  coupling  between  Eqs.  (5)  and  (6)  is  through  the  viscosity  ratio, 
N,  which  is  a strong  function  of  temperature. 

Since  we  are  interested  in  the  conditions  for  incipient  separa- 
tion, we  seek  solutions  to  the  above  set  that  also  correspond  to 
f" (0)  * 0.  These  solutions  would  then  show  the  vanishing  of  the  shear 
stress  at  the  wall,  which  is  the  usual  definition  of  separation. 

If  the  Prandtl  number  is  large,  then  the  thermal  boundary  layer 
is  much  thinner  than  the  velocity  boundary  layer.  This  suggests  there 
is  a thin  region  near  the  wall  where  vircosity  variation  is  important 
in  the  momentum  balance.  For  such  a layer,  the  dominant  terms  in  the 
force  balance  are  shear  and  pressure  gradient;  inertial  effects  are 


negligib le. 

We  introduce  a new  variable,  Z,  such  that  n = eZ,  so  Z * 0(1) 
when  n = 0(e);  the  dependence  of  this  small  parameter,  e,  on  physical 
quantities  is  determined  later.  Introduce  a new  stream  function  vari- 
able, X(Z),  such  that 


fnnn(n)  * exzzz(z)  * 

(7) 

fnn(n)  = e6xzz(z>  * 

(8) 

f (n)  = e2bx_(z)  , 
n z 

(9) 

f (n)  = e3Bx(z)  , 

(10) 

(0)  = f'(0)  = f(0)  have  been  used. 

This  order- 

ing  of  f in  terms  of  e follows  from  the  assumption  that  N is  a quantity 
of  0(1). 

Then  various  terms  in  the  momentum  equation  are  of  the  order: 


« Vn  + f£nn  + 8(1  ' V * 0 

0(1)  0(e4)  0(1  - e4) 


The  determination  of  e follows  from  the  consideration  of  the  energy 
equation;  where  g is  of  0(1) 


g + Pr  f g =0 

00  n 


°fe)  “M6) 


It  if  necessary  that  there  be  a balance  between  conduction  and  trans- 
port if  the  boundary  conditions  are  to  be  satisfied  at  the  wall  and 
the  edge  of  the  thermal  layer.  From  this  balance,  the  order  of  magni- 
tude  of  e is  established:  e = (-B  Pr^)  . The  minus  sign  before  8 is 
chosen  because  B is  a negative  quantity.  While  this  quantity,  e,  is  not 
particularly  small  if  Pr^  ~ 10  and  -B  ~0.2,  there  is  empirical  evi- 
dence from  the  work  of  Liepmann^  and  Narasimha  and  Vasantha^  that 
even  if  Pr  is  as  small  as  0.7,  solutions  obtained  via  the  approxima- 

OO 

tion  e ■*  0 are  very  accurate.  The  reduced  momentum  equation  becomes 


(NXZZ>Z  + 1 ■ 0 


(13) 
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and  the  reduced  energy  equation  becomes 

gZZ  " XgZ  = 0 • 


(14) 


In  general,  these  two  nonlinear  equations  must  be  solved  simultaneously; 
they  are  coupled  through  the  dependence  of  the  viscosity  ratio,  N,  on 
g,  and  the  dependence  of  g on  X.  However,  this  solution  would  be  easy 
to  obtain  numerically.  The  appropriate  boundary  conditions  on  X are 
that  X(0)  = Xz (0)  = 0,  and  the  conditions  on  g are  that  g(0)  = 1,  and 
g('/  - 0.  The  condition  that  X (0)  = 0 follows  from  the  additional 

requirement  that  the  relationship  between  B and  T - T be  such  that 

w 00 

the  shear  vanishes  at  the  boundary. 

While  it  would  be  interesting  to  solve  the  coupled  nonlinear  set 
(Eqs.  (13)  and  (14))  numerically,  we  can  proceed  to  obtain  quantitative 
estimates  by  introducing  an  additional  approximation  in  the  energy 
equation.  Prior  to  stating  that  approximation,  however,  it  is  necessary 
to  integrate  the  momentum  equation. 

The  first  integral  of  the  momentum  equation  yields 


NX 


ZZ 


-z 


» 


the  next  integral  yields 


X 


Z 


and  the  last  integral  yields 

Z 


Z dZ 
N (Z) 


Z”  dZ" 
N(Z”)  • 


(15) 


(16) 


X 


(17) 
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We  integrate  Eq.  (16)  by  parts 


In  the  limit  as  Z •+■ 


CD 


where  we  anticipate  that  N + 1 as  Z + 


(18) 


(19) 


VISCOSITY  VARIATION  WITH  TEMPERATURE 

It  is  now  necessary  to  specify  N,  the  viscosity  ratio,  as  a func 
tion  of  g.  Gazley*  has  shown  that  a good  approximation  for  the  kine- 
matic viscosity  of  water  in  the  range  of  temperature  between  40°F  and 
100°F  is 


V 


10~5 

.0807  + 0.0126  T 


and  between  100°F  and  160°F  it  is 


(20) 


V 


10^ 

-0.1946  + 0.0153  T 


(21) 


2 

where  T is  in  °F  and  v is  in  ft  /sec.  This  suggests  a viscosity- 
temperature  model  of  the  form 


1 

a + bT  ' 


(22) 


*Personal  communication  from  Carl  Gazley,  Jr.,  The  Rand  Corporation, 


1 


1975. 
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r 


For  such  a model. 


I a + bT 
N ’ a + bT, 


a 

a + bTm 
OO 


+ 


b (X  + <TW  - TJs] 

* + bTw 


= 1 + oATg 


(23) 


! 


where  AT  = T - T^,  and  a = b/ (a  + bTj  . 

w -1 
For  water  in  the  temperature  range  between  40°F  to  100°F,  a « 0.153  (°F) 

Using  the  relationship 


the  reduced  velocity  function  dX/dZ  is  written 


dX 

dZ 


dZ 


(24) 


(25) 


APPROXIMATE  SOLUTION  OF  THE  ENERGY  EQUATION 

It  is  traditional  in  heat  transfer  analyses  to  employ  crude  approx- 
imations to  the  convective  terms  in  the  energy  equation  and  still  obtain 
highly  accurate  results.  Our  approach  continues  that  tradition,  and  we 
trust  that  further  precise  solution  of  Eqs.  (13)  and  (14)  will  confirm 
our  approximate  analysis. 

We  assume  that  for  purposes  of  solving  the  energy  equation  a suit- 
able  approximation  to  X is  X = -Z  /6N,  where  N is  e average  reference 
value  of  N evaluated  between  T^  and  T . (See  Appendix  A for  a more 
vigorous  justification  of  this  approximation.)  This  is  equivalent  to 
using  the  constant  property  version  of  X,  where  the  viscosity  is  evalu- 
ated at  some  reference  condition  in  the  boundary  layer.  Notice  that 


I 

j 


-Z 

x ■*  - — at  large  values  of  Z . 


(26b) 


This  is  equivalent  to  approximating  f(n)  by 


f (n)  - - 

3N 


(27) 


in  the  original  variables. 

The  resulting  energy  equation  is 


4 + . 

dZ  6N 


(28) 


whose  solution  is 


L 

/. 


E = A I exp(-Z  /24N)  dZ  + C . 

'o 


(29) 


The  boundary  conditions,  g(0)  * 1,  g(°°)  * 0,  fix  the  two  constants  A 
and  C so  that 


g(Z)  = 1 - 


L 

/ 


exp(-Z  /24N)  dZ 


f. 


exp(-Z^/24N)  dZ 


(30) 


I 

t 
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Therefore , 


d_ 

dZ 


uAT 


OO 


• exp(-Z^/24N) 
exp(-Z4/24N)  dZ 


(31) 


and  the  asymptotic  expression  for  dX/dZ  (Eq.  (19))  can  now  be  evaluated. 
As  Z + the  velocity  function 


2 

d X ^ Z 
dZ  ‘ 2 


aAT 

2 


Z2  exp(-Z4/24N)  dZ 


exp(-Z4/24N)  dZ 


(See  Appendix  B for  evaluation  of  integrals.) 
Thus 


(32) 


or 


aAT  <a^r(3/4)| 

- ru/4)} 


2 

d x K zr 

dZ  2 


aAT  f(3/4) 

2 r(i/4) 


(24N)1/2 


or 


dX 

dZ 


Z 

2 


(24N)"1/2  c 


(33) 


(34) 


(35) 
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\ 

1 


where  c - F(3/4)/T(l/4) . It  is  important  to  note  that  thr  constancy 

of  the  second  term  on  the  right-hand  side  of  Eq.  (35)  is  not  a conse- 

3 

quence  of  the  approximation  X = - Z /6N,  and  that  because  of  the  in- 
equality, Eq.  (A. 8)  of  Appendix  A,  the  second  term  on  the  right-hand 
side  of  Eq.  (35)  must  be  a constant  with  lower  and  upper  bounds  given 
by 

CD  2 

~ (24Nw)'1/2  c f Z2  exp  f X(Z')  dZ'j  dZ  < ~ (24)~1/2  c 

*0  ' •'o  ' 


(see  Eq.  (A. 3)). 

APPROXIMATE  SOLUTION  OF  THE  OUTER  LAYER 

We  have  thus  determined  the  structure  of  the  velocity  profile 
within  the  thermal  boundary  layer.  We  recast  the  asymptotic  formula 
for  dX/dZ  in  terms  of  the  original  variables,  n and  f: 


df 

dn 


Be 


2 dX 
dZ 


(Z) 


(36) 


= B(-S  § (Z) 


(37) 


or 


or 


df 

dn 


B(-B  PrJ 


-1/2 


[-* 


otATc 

2 


(24N)1/2 


(38) 


^ = - f n2  - f (-B  PrJ-172  aATc(24N)1/2 


(39) 


Thus,  as  the  variable  n 0,  the  outer  flow  must  be  of  the  form 


u 

U 


- f n2  + k , 
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where  ^ is  a positive  constant,  in  order  to  match  to  the  velocity 
profile  at  the  outer  edge  of  the  thermal  boundary  layer. 

Stewartson(8)  has  examined  solutions  of  the  Falkner-Skan  equations 
that  have  this  property.  They  are  two-dimensional  wake-like  flows  that 
have  a positive  slip  velocity  at  n - 0 as  well  as  zero  shear  and  possess 
values  of  8 between  -0.198  and  -0.5.  Berger (9)  discusses  these  flows 
and  graphs  the  relationship  between  df/dn  (0)  and  8.*  Thus  it  is  only 
necessary  to  use  that  relationship  between  ^(0)  and  8,  in  addition  to 
e9-  (39)  evaluated  as  n -*•  0, 


5H  ’ -6C-6  rr „>l«  (24S)1'2 


(40) 


in  order  to  determine  the  relationship  between  AT  and  8 that  corre- 
sponds to  incipient  separation.  (See  Fig.  1 for  a sketch  of  the  outer 
and  inner  velocity  profiles.) 

Figure  26  of  Ref.  9 is  replotted  in  the  convenient  form 


df 


dn 


(0)/|B|1/2  versus  8 


This  is  shown  in  Fig.  2. 
Since 


| »)/|6|1/2  - <Pr_)-l«  2£c  (245)l/2  i 


the  ordinate  is  proportional  to  (T  - T ) , 


NUMERICAL  VALUES 
At  T - 60*F, 


- 1.22  x lo"5  ft2/sec  , 


b - 0.0126  x 105(ft2/sec)_1(°F)_1 


See  Fig.  26,  p.  75,  of  Ref.  9. 


hii 


■ 


(41) 


NOTE:  The  outer  velocity  profile,  if  continued  to  the  wall,  would 
follow  the  dashed  line  and  would  correspond  to  a slip  velocity, 
Kj,  and  zero  shear.  Thn  Inner  velocity  profile  corrects  this  and 
brings  the  velocity  at  y=0  to  zero  while  still  maintaining  shear 
at  y=  0.  The  dashed  line  is  a parabola  (see  Eq.  (3^)  and 
corresponds  to  the  effect  of  the  surface  heating  on  the  outer 
velocity  profile. 

Fig.l  — Sketch  of  the  velocity  profile 

and 

a = 0.0153  (°F)-1  , 

where  the  density  variation  is  neglected. 

If  Pr  w 8,  then 

(Pr  )"1/2  ?£  (2AN)1/2  AT  = 0.00453  AT(N)1/2  . ( 


F"  (o) 


a ATc 


NOTE:  7-  m/5  (24  N) 


w* 


1/2 


= 0.0045  AT  if  Too  = 60°F 
and  N = 1 


Fig.  2 — Slip  velocity  versus  Falkner-Skan  P 
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111.  RESULTS  AND  CONCLUSIONS 

The  required  wall  overheat,  AT  . , is  the  minimum  value  of  the 

min 

temperature  difference,  T^  - T^,  that  will  prevent  separation  for  a 
particular  negative  value  of  0.  Its  dependence  on  0 is  shown  in  Fig. 
3,  where  60°F  water  is  the  ambient  fluid.  This  figure  is  based  on 
assuming  a reference  value,  N,  of  the  viscosity  ratio  N = y/yro>  which 
corresponds  to  unity.  Since  N decreases  with  increasing  temperature, 
any  other  choice  for  the  reference  value  N would  be  less  than  1.  From 
Eq.  (42),  it  is  easy  to  include  any  other  choices  for  N,  such  as 


which  corresponds  to  an  average  reference  viscosity  ratio.  The  straight- 
line  relationship  between  ATmin  and  0 is 

AT  . « -1000  {0  + 0.198]  . (43) 

mm 

When  an  average  reference  condition  is  used,  this  straight  line  becomes 
one  with  positive  curvature,  and  the  resulting  value  of  AT  ^ is  always 
greater  than  the  AT^  which  is  computed  for  N equal  to  one.  For  ex- 
ample, if  ATmin  is  40°F,  Fig.  3 indicates  that  0 = -0.24,  and  the  choice 
of 


leads  to  a value  of  0 which  is  about  -0.235.  This  results  because  N 
is  0.9  instead  of  unity. 

From  Fig.  3,  we  see  that  a wall  overheat  of  40°F  yields  a pressure 
gradient  parameter  0 which  is  about  -0.24  at  separation.  This  corresponds 
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Fig.  3 — Required  surface  heating  and  suction  to  delay 
laminar  separation  beyond  p=  -0.198 
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to  a 25  percent  increase  in  magnitude  over  the  classic  ^alkner-Skan 
value  of  B = -0.198.  It  is  easy  to  relate  this  to  the  pressure  grad- 
ient or  velocity  gradient  using  the  following  definition  of  B: 


If  separation  occurs  aft  of  the  maximum  diameter,  then  U ~U  . r can 

0 00  * 

be  approximated  as  a constant,  and  the  major  variation  in  B comes  from 
the  velocity  gradient  dUg/dx.  The  fractional  increase  in  the  magnitude 
of  B then  corresponds  roughly  to  the  fractional  increase  in  the  magni- 
tude of  dU  /dx. 
e 
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Thus,  a 25  percent  increase  in  the  value  of  8 corresponds  to  a 

similar  increase  in  the  permitted  magnitude  of  the  velocity  gri  ient 

without  separation.  For  the  class  of  exact  wedge  flows,  the  free- 

stream  velocity  is  now  U ~x‘0,11  rather  than  U ~ x"0.091  as  it  is 

e e 

for  the  constant  property  case. 

In  order  to  illustrate  the  effectiveness  of  suction  compared  to 
wall  heating  in  delaying  separation,  we  have  shown  in  Fig.  3 the  re- 
quired suction  velocity  ratio,  v/U^,  which  will  prevent  separation  for 
the  same  class  of  Falkner-Skan  'lows.  The  quantity  U x/v  is  the 
Reynolds  number  and 

J uer2;1  dx 
A H * 

e 

For  a boundary  layer  near  the  maximum  diameter  point,  the  quantity  A 
is  close  to  unity. 

The  suction  velocity  ratio,  v /U  , which  results  in  the  same  ef- 

w 00 

feet  on  separation  as  a 40°F  wall  overheat  is  about 


0.1 


If  the  arc  length  Reynolds  number  at  the  maximum  diameter  is  10  x 106, 

_ O 

then  the  suction  velocity  ratio  is  0.03  x 10  ; i f the  Reynolds  number 

6 o 

at  the  maximum  diameter  is  20  x 10  , the  value  of  v /U  is  0.02  x 10  . 

w °° 

From  Fig.  3,  it  is  clear  that  a small  amount  of  suction  has  a very 
powerful  effect  on  the  delay  of  separation.  Figure  4 is  taken  from 
Rosenhead^^  and  shows  that  enormous  unfavorable  pressure  gradients 
(8  < -1)  can  be  sustained  without  separation.*  This,  of  course,  is  in 
accord  with  both  intuition  and  experiment.  Wall  heating,  as  we  have 
shown,  has  a much  less  dramatic  effect,  and  it  requires  wall  overheats, 

* 

See  Fig.  V.10,  p.  249,  of  Ref.  10. 


>1=3 


/Uer2jdx_ 

U#r,J* 

j = 0,  2-Dim 
I , A <i  -tym 


Folkner-Skon  pressure  gradient  parameter,  p 


Fig.  4 — Suction  required  to  delay  laminar  separation  beyond  p=  -0.198 


Tw  _ T»»  whlch  are  unrealistically  large,  to  have  such  significant 
effects. 

The  analysis  we  have  shown  here  is  an  approximate  one,  but  the 
approximations  are  well  understood  and  are  not  expected  to  degrade  the 
numerical  results  significantly.  However,  we  have  examined  in  detail 
only  the  Falkner-Skan  separating  flows.  Both  approximate  and  precise 
numerical  results  for  other  more  realistic  flows  should  be  obtained, 
since  the  weakest  part  of  our  approach  is  the  conceptual  extension  of 
our  results,  using  a crude  local  similarity  argument,  to  non-Falkner- 
Skan  flows. 

Our  analysis  predicts  that  laminar  separation  can  be  delayed  sig- 
nificantly if  the  adverse  pressure  gradient  is  small  in  the  region  of 
laminar  separation.  Nevertheless,  it  appears  that  surface  heating  in 


-20- 


water  has  much  less  Impact  on  delaying  separation  than  does  even  a small 
amount  of  suction. 


f 


■l 
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Appendix  A 

the  APPROXIMATION  X ~ -z3 


X ~ - 


A^ore  rigorous  justification  of  the  use  of  the  approximation 
Z in  the  energy  equation  can  be  given  as  follows. 


The  matching  rule  used  in  Eqs.  (36)  through  (39)  is  formally 


lim  df (n) 
TY*0  dn 


1 im  d , 3 

z-*»  dn  e 6X(Z)1 


2 . lim  dx 
Z-x»  dZ 


or 


torn' 

\ dn  , 


B (i) 

' 'a 


(A.  1) 


where  (dX/dZ)^  is  determined  from  Eq.  (25).  Equation  (25),  which  is 
exact,  requires  that  dg/dZ  be  known  from  the  solution  of  the  energy 
equation.  In  this  connection,  we  note  that  the  energy  equation,  Eq. 
(14),  can  formally  be  solved  exactly  for  arbitrary  X(Z), 


d& 

dZ 


-B 


exp 


/ Z 

(/  x(z)  dz) 


(A.  2) 


where  B is  an  arbitrary  real  constant  determined  by  the  boundary  con- 


ditions on  g(Z),  and  the  constant  in  Eq.  (A. 2)  is  written  as  B2  to 


emphasize  the  fact  that  dg/dZ  < 0 for  all  Z. 
into  Eq.  (25),  we  obtain 


Substituting  this  result 


dX 

dZ 


zr_ 

2N 


aB2AT 


fQ  22  CXP  a X(Z')dZ' j 


dZ 


(A.  3) 
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It  is  the  behavior  of  this  expression  as  Z +“  that  we  are  interested 
in. 

Consider  the  integral  appearing  in  Eq.  (A. 3).  We  must  first  begin 
by  looking  at  the  expression  for  dg/dZ  appearing  therein.  We  note  that 
it  follows  from  the  fact  that  f(Z)  >0  for  all  Z and  the  definition  of 
X(Z)  that  X(Z)  < 0 for  all  Z.  It  then  follows  for  any  B that  dg/dZ 
approaches  zero  as  Z -*■  °°.  We  then  have  that  dg/dZ  is  an  exponentially 
decaying  function  of  Z for  large  Z and  that  g(Z)  decreases  monotonically 
from  1 at  the  wall  to  zero  at  infinity. 

For  the  purposes  of  this  report  it  is  sufficient  that  we  have  bounds 
on  the  second  term  on  the  right-hand  side  of  Eq.  (A. 3).  This  can  be 
done  by  analyzing  Eq.  (A. 2)  more  carefully.  First,  recall  that 


X(Z)  ~- 


near  Z * 0 , 


(26a) 


X(Z)  ~ - r 


for  Z -*•  00  , 


(26b) 


and  that  with  the  assumption  that  AT  * Tw  - > 0,  it  follows  from  the 

definition  of  N (see  Eq.  (23))  that  Nw  < 1.  Thus  -Z  /6Nw  < -Z  /6. 
Further,  from  Eq.  (16)  we  see  that 


Xz  < 0 


for  all  Z 


(A. 4) 


Also,  since  from  Eq.  (23) 


1 dN  d£ 

- ^ dZ  * 041  Jz 


it  follows  from  dg/dZ  < 0 that  dN/dZ  > 0.  This  in  turn  means  that  N(Z) 
is  monotonically  increasing,  and  in  particular  it  satisfies 


N s:  N (Z)  £ 1 

W 


(A. 5) 
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If  we  use  this  last  inequality,  it  follows  immediately  from  Eq. 
(16)  that 


-r/ f 

w 0 */o  •'o 


or 


2N 

w 


s:  X„ 


(A. 6) 


We  also  note  that  since  N - > 0,  it  follows  from  Eq.  (15) 

that 


Xzz  < 0 .or  all  Z . (A. 7) 

The  consequence  of  Eqs.  (A. 4),  (A. 6),  and  (A. 7)  is  that  the  curve 

3 3 

X(Z)  must  lie  between  the  two  curves  -Z  /6N  and  -Z  /6,  as  shown  in 

w 

Fig.  A.  1,  coinciding  with  the  first  of  these  near  the  origin  and 

approaching  the  latter  as  Z -*  °°.  For,  since  X(Z)  starts  at  the  origin 

with  the  same  value,  slope,  and  curvature  as  -Z  /6N  , it  can  "get  out- 

w 

side"  this  curve  only  if  it  violates  the  left-hand  part  of  inequality 
Eq.  (A. 6)  or  doubles  back  on  itself.  The  latter  of  these  alternatives 
is  also  ruled  out  since  it  would  lead  to  X(Z)  being  a double-valued 

'i 

function  of  Z.  On  the  other  hand,  X(Z)  cannot  "get  outside"  -Z  /6,  for 
to  do  so  and  still  satisfy  the  right-hand  part  of  inequality  Eq . (A. 6) 
would  require  that  it  double  back  on  itself  and  again  be  double-valued. 
Thus  X(Z)  is  bounded  as  shown  in  Fig.  A.l. 

From  Fig.  A.l  we  note  that 


(A. 8) 


i 


where  N is  some  average  reference  value  of  N evaluated  for  some  T lying 
between  and  T . This  is  the  approximation  made  in  the  main  text. 
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Hovever,  it  should  be  noted  that  for  the  purposes  of  this  report  what 
is  required  is  the  minimum  value  of  the  temperature  difference  T - T 

* W oo* 

which  will  have  a desired  effect,  the  magnitude  of  which  effect  is 
determined  by  the  value  of  the  right-hand  side  of  Eq.  (A. 3)  for  that  AT. 
This  is  obtained  from  Eq.  (A. 8)  by  assigning  to 


XCZ'JdZ1 


its  upper  limit,  namely  -Z  /24.  In  fact,  the  specific  numerical  esti- 
mates given  in  the  report  involve  the  use  of  this  latter  value,  since 
they  are  based  on  N « 1, 
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Api end lx  B 


EVALUATION  OF  INTEGRALS 


Two  integrals,  which  must  be  evaluated  (see  Eq.  (32)),  are 


I : / Z2  exp(-ZA/24N)  dZ 

** r\ 


I 2 * / exp(-Z/24N)  dZ  . 


We  consider  1^  in  detail.  Let 


(B.  1) 


Z = (24N)1/4  t1/A  , 


(B.  2) 


d2  . t-3/«  at  . 

4 


(B.  3) 


Thus,  we  can  write  1^  as 


(24N)3/4 


t-!/4  e-t  dt 
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Slnce 


» 


then 


Similarly, 


so 


r(n) 


00 


■■  ro/4) 


r(i/4)  . 


1 

I 

(B.4) 


t dt  (B.5) 


(B.6) 


l 

j 
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